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Q ' I describe the Einstein's gravitation of 3+1 dimensional spacetimes using the (2,2) formal- 

* ^ • ism without assuming isometries. In this formalism, quasi-local energy, linear momentum, 

I ■ and angular momentum are identified from the four Einstein's equations of the divergence- 

type, and are expressed geometrically in terms of the area of a two-surface and a pair of null 
^►^ I vector fields on that surface. The associated quasi-local balance equations are spelled out, 

f^ — ' and the corresponding fluxes are found to assume the canonical form of energy-momentum 

^ ' fiux as in standard field theories. The remaining non-divergence-type Einstein's equations 

\^^ . turn out to be the Hamilton's equations of motion, which are derivable from the non- 

^P , vanishing Hamiltonian by the variational principle. The Hamilton's equations are the evo- 

f— ^ I lution equations along the out-going null geodesic whose ajfine parameter serves as the time 

^^ ' function. In the asymptotic region of asymptotically flat spacetimes, it is shown that the 

' quasi-local quantities reduce to the Bondi energy, linear momentum, and angular momentum, 

5-H i and the corresponding fluxes become the Bondi fluxes. The quasi-local angular momentum 

. . , turns out to be zero for any two-surface in the flat Minkowski spacetime. I also present a 

. !^ I candidate for quasi-local rotational energy which agrees with the Carter's constant in the 

^\f ' asymptotic region of the Kerr spacetime. Finally, a simple relation between energy-flux and 

j^ ■ angular momentum-flux of a generic gravitational radiation is discussed, whose existence 

reflects the fact that energy-flux always accompanies angular momentum-flux unless the flux 
is an s-wave. 
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I. INTRODUCTION 

It has been known for sometime that the Plebahski equation|]J, the self-dual Einstein's equation 
of 4 dimensional Euclidean space, can be obtained as the large n limit of the equations of motion 
of a certain class of s/(n)-valued non-linear sigma models in 2 dimensions [3, y,|j|- The equivalence 
of these equations defined in two different dimensions is quite unexpected, but if one realizes that 
large n limit of the sl{n) Lie algebra is iust the Lie algebra of area-preserving diffeomorphisms 
of an auxiliary 2 dimensional surfacejsl, |g, |3], and that the equations of motion of s^(oo)-valued 



* E-mail address :yoonjh@konkuk.ac.kr 



non-linear signia models in 2 dimensions are in fact partial differential equations on 4 dimensional 
space, then one might be more comfortable with the idea of describing 4 dimensional self-dual 
Einstein's gravity as a limit of a certain class of 2 dimensional field theories, and can show that 
the two theories are in fact identical. This correspondence is supported further by the observation 
that the Plebahski equation and the s/(n)-valued 2-dimensional non-linear sigma models are both 
integrable. 

One may be interested in extending this idea of describing 3+1 dimensional theories from 1+1 
dimensional perspective without the self-dual restriction and in the Lorentzian regime. There are 
several advantages of such a description, if it is possible at all. They stem from the fact that 
1+1 dimensional field theories are usually more manageable than 3+1 dimensional ones, both 
classically and quantum mechanically. For example, a number of field theories in 1+1 dimensions 
are renormalizable, due to the dimensionlessness of field variables in a naive power counting. There 
would be an enormous gain if one ever succeeds in describing the Einstein's gravitation in 3+1 
dimensions as a limit of some kind of 1+1 dimensional field theories which eventually proves to be 
renormalizable. This idea sounds strange but does not seem impossible, since the renormalizability 
is highly sensitive to the spacetime dimensions on which the theories are defined. 

These reasonings led us to seek the possibility whether the Einstein's gravity in 3+1 dimensions 
without the self-dual restriction is describable as a 1+1 dimensional field theory [3]. The idea was 
simply to split a 3+1 dimensional spacetime into a 1+1 dimensional base manifold and a 2 dimen- 
sional fibre space, and write down the Einstein-Hilbert action. Then the Einstein-Hilbert action 
becomes 1+1 dimensional field theory action, where the infinite dimensional group of diffeomor- 
phisms of 2 dimensional fibre space becomes the Yang-Mills gauge symmetry. But this program 
was successful only in a formal sense, since the resulting 1+1 dimensional action did not seem much 
useful, which made the whole idea of describing the Einstein's gravitation as 1+1 dimensional field 
theory questionable. The follow-up idea was to use the gauge freedom of the 3+1 dimensional 
spacetime y, |lO|, |ll|, |l^, |lj, |lj|. If one chooses one of the spacetime coordinates as the affine 
parameter of the out-going null geodesic, then it turns out that the 1+1 dimensional field theory 
description of the Einstein's gravitation is simplified significantly. 

The purpose of this paper is to present several unexpected results that I obtained in the (2,2) 
fibre bundle description of the Einstein's gravitation, and discuss their physical implications. First, 
I will present quasi-local balance equations of energy, linear momentum, and angular momentum 
for an arbitrary compact two-surface, which are just two-surface integrals of the four divergence- 
type equations that are part of the Einstein's equationsjlj, |lj, |l5|. Quasi-local energy, linear 
momentum, and angular momentum are expressed in the coordinate-independent and geometric 
way in terms of the area of a two-surface and the in- and out-going null vector fields at each point 
of that surface lalldl- They are Bondi-like, since their rates of changes are given by fiuxes of the 
canonical form lq| 

roa??"~^^/£y, (1.1) 



where t] is an appropriate vector field defined at each point of a two-surface. 

Second, problems of defining quasi-local an gula r momentum and associated rotational energy 
have been particularly subtle issues [l^l, |2IJ) UM, l2^; Uai l2J|- This is due to the fact that the very 



notion of rotation depends on the choice of the coordinates, which implies that one can always 
remove the effects of rotation by working in a co-rotating coordinate system. On the other hand, 
it is natural to demand that the angular momentum and the rotational energy of any compact 
two-surface in the flat Minkowski spacetime be zero. It will be seen that our quasi-local angular 
momentum and rotational energy not only become zero for any two-surface in the flat Minkowski 
spacetime, but also reduce to the standard values of the total angular momentum and the Carter's 



constant in the asymptotic region of the Kerr spacetime, resDectivelv|25l. l26l. l27l. l28l| . In this sense, 
our quasi-local rotational energy may be regarded as a quasi-local generalization of the Carter's 
constant of a generic gravitational field. 

Third, using the affine parameter of the out-going null geodesic as the time coordinate, I will 
write down the Hamiltonian of the Einstein's theory [22|. I will obtain the Hamilton's equations of 
motion from this Hamiltonian using appropriate boundary conditions, which determine the time 
flows of the field variables. Together with the quasi-local balance equations (or the constraint 
equations depending on the signature of the 3-dimensional hypersurface) , it will be seen that the 
Hamilton's equations of motion constitute the full Einstein's equations. 

Finally, I will present a simple but general relation between quasi-local energy-flux and angular 
momentum-flux of a generic gravitational radiation that has no isometries. It is a generalization 
of the well-known relation of mass-loss and angular momentum- loss |3Cll| . 

6U = —6L^ (1.2) 

rriz 

for small perturbations around a stationary and axi-symmetric spacetime, where u) and m^ are the 
frequency and azimuthal angular momentum of the perturbations, respectively. To my knowledge, 
such a relation between these gravitational fluxes of the most general type has not been discussed 
before, but it strongly indicates that our identifications of fluxes are physically correct, since 
energy-flux always carries angular momentum-flux unless the radiation is an s-wave. 

This paper is organized as follows. In section^ I will introduce the kinematics of the (2,2) fibre 
bundle formalism, and write down the Einstein's equations. Then I will discuss 1+1 dimensional 
gauge theory aspects of the Einstein's gravitation of 3+1 dimensions from this fibre bundle point 
of view. 

In section 11111 I will study the four Einstein's equations that are first-order in the derivatives 
along the out-going null vector field. These equations, which are the natural analogs of the Ein- 
stein's constraint equations in the 3+1 formalism, turn out to be divergence- type equations. It is 
from the two-surface integrals of these equations that one obtains quasi-local balance equations 
of gravitational energy, linear momentum, and angular momentum. I will also present quasi-local 
gravitational rotational energy. The Carter's constant, which is usually interpreted as a measure of 
intrinsic rotation of gravitational field, is known to exist for a certain class of spacetimes that have 



two commuting Killing symmetries, and for the Kerr spacetime, it is just the total angular momen- 
tum squared. Our quasi-local rotational energy reduces to the Carter's constant for asymptotically 
Kerr spacetimes, as is shown in IVII and therefore, may be regarded as a quasi-local generalization 
of the Carter's constant to spacetimes that have no isometries. 

In section IIVI it will be shown that the remaining Einstein's equations, which are second- 
order in the derivatives along the out-going null vector field, are the Hamilton's equations of 
motion derivable from a non-vanishing Hamiltonian by the variational principle. The details of 
this derivation are given in Appendix. Thus, together with the quasi-local balance equations (or 
constraint equations depending on the signature of the 3-dimensional hypersurface) , the Hamilton's 
equations of motion make up for the full Einstein's equations in this formalism. 

In section El quasi-local energy, linear momentum, and angular momentum of the previous 
sections will be expressed in the coordinate-independent and geometric way, using the area of 
a two-surface and a pair of null vector fields orthogonal to that surface. Relative to a given 
background spacetime against which these quasi-local quantities are measured, quasi-local energy 
and linear momentum are given by the rates of changes of the area of the two-surface along the in- 
and out-going null vector fields, respectively, and quasi-local angular momentum associated with a 
vector field S, is given by two-surface integral of the projection of the twist of the in- and out-going 
null vector fields onto ^ modulo a background-dependent subtraction term. 

In section IVII I will study the quasi-local balance equations at the null infinity and show that 
they all agree with the well-known Bondi formulae of energy-loss, momentum- loss, and angular 
momentum-loss. In order to show these correspondences, it is necessary to find the asymptotic 
fall-off rates of the metric and their derivatives near the null infinity, using the affine parameter 
of the out-going null geodesic as the radial coordinate [ij, Il4l . l3ll . l32l . l3.3i | . I will present the 
asymptotic fall-off rates in this section. It will be shown that the quasi-local rotational energy in 
the asymptotic region of the asymptotically Kerr spacetimes agrees with the Carter's constant of 
the Kerr spacetime. 

In section rV 111 a general relation between quasi-local energy-fiux and angular momentum-flux 
will be presented for a generic gravitational radiation. When restricted to small perturbations 
around a stationary and axi-symmetric spacetime, it will be shown that this relation reduces to 
the well-known relation of mass-loss and angular momentum-loss in the perturbation theory of the 
Kerr black hole 



3C|. 



In the Appendix, I present in detail the derivation of the non-divergence type Einstein's equa- 
tions as the Hamilton's equations of motion associated with a non-vanishing gravitational Hamil- 
tonian. 

II. KINEMATICS 

da = d/dy"-, d+ = d/du, (9_ = d/dv 

In this section, I will introduce the kinematics of the (2,2) fibre bundle formalism 
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and 





d+ := d+ 


- A^da, 




5_ := a_ 


- A^da, 


d 


^- = IT' 
at" 


dy' 



write down the Einstein's equations. This section serves mainly to fix the notations. Let us consider 
the following line element 

ds^ = -2dudv - 2hdu^ + (t)ab {dy" + Aj^du + A^dv) (dy^ + A^du + A^dv) , (2.1) 

where +, — stands for u,v, respectively J9|, |]^, |l]J, |lj, UJ, |lj, |3a, |37|, |3a, |39|] . In order to understand 
the geometry of this metric, it is convenient to introduce the following vector fields {d±} defined 
as 

(2.2) 
(2.3) 

where 

The inner products of the vector fields {d±,da} are given by 

< d+, d+ >= -2h, < d+, d- >= -1, < d-, d- >= 0, 

< 4, da >= 0, < da, db >= CPab- (2.5) 

The hypersurface u = constant is an out-going null hypersurface generated by d- whose norm is 
zero. The hypersurface v = constant is generated by 3+ whose norm is — 2/i, which can be either 
negative, zero, or positive, depending on whether 5+ is timelike, null, or spacelike, respectively. 
The vector fields {d±} are called horizontal since they are orthogonal to {da}, and two dimensional 
section spanned by {d±} has the Lorentzian signature. The intersection of two hypersurfaces 
u,v = constant defines a spacelike two-surface A'^2 labeled by {y'^}, which is assumed to be compact 
with a positive-definite metric 4>ab on it (see FIG. 1). The metric (pat is decomposed into the area 
element e"^ and the conformal two-metric pab normalized to have a unit determinant 

(l^ab = e'^Pab (det Pab = 1). (2.6) 

For later uses, let us express the in-going null vector field n and out-going null vector field / in 
term of {d±}. They are given by 

n = d+-hd^, (2.7) 

/ = a_, (2.8) 

and satisfy the normalization condition 

< n, l>= -1. (2.9) 

Notice that d/dv is either spacelike or null, since its norm is given by 

< I;, I; >= e^abA.^A^ > 0. (2.10) 



l = dldif(a = 2,i) 



2 dimensional 
fibre space A^2 




1+1 dimensional 
base manifold 



FIG. 1: This figure shows the geometry of the (2,2) fibre bundle sphtting of 3+1 dimensional spacetime. 
The 1+1 dimensional base manifold is spanned by {d±} and the 2 dimensional fibre space N2 by {da}- The 
horizontal vector fields {d±} are orthogonal to 7V2, and A^ are the connections valued in the Lie algebra of 
the difl^eomorphisms of N2 ■ 



The coordinate v increases uniformly as I evolves, since we have 

£iv = 1. 
In the gauge where A_^ = 0, / is given by 

I 



d_ 

dv'' 



(2.11) 



(2.12) 



which tells us that v becomes the afRne parameter of the out-going null geodesic I. 



The complete set of the vacuum Einstein's equations are found to be 



S){ 



(a) e'^D+D^a + e^D^D+a + le" {D+a){D^a) - le" {D^h){D^a) - -e^" pabF^^F^" 



+e-i?2 - he''{{D^af - ip'^V''(i^-Pac)(i^-P6d)} = 0, 



(2.13) 



(b) 



-e^Dla - -e^iD+af - e''{D^h){D+a) + e" {D+h){D^(j) + 2he''{D^h){D^a) 

+h{ - e^D+a){D^a) + ^e" p-'p'\D+pa,){D_pM) + ^e^V^.F+^F^^ - e^R^} 
+h\'^{{D^a)^ - i/V'(^-Pac)(I)-P6d)} = 0, (2.14) 

(c) 2e''{Dla) + e"{D^af + ^e>"V''(^-/0ac)(I?-P6d) = 0, (2.15) 



= 0, (2.16) 

(e) -I)+(e2>,bF+^) - e''da{D+a) - ^e>V(^+P6d)(5aPce) + db{e'' p'^'D+pac 

+2he''da{D_a) + he^p'"'p''%D_pM){daPce) + 2e'^5„(D„/i) - 2db[he'' p''D_pac 
= 0, (2.17) 

(/) -2e''D'ih - 2e''{D_h){D_a) + e^D+D^a + e^'D^D+a + e''{D+a){D^a) 

+ le^p-^p'^\D^Pac){D^PM) + e2>,,F+_"F+^ - 2h^^ [ola + ^(I^-ct)^ 

+\p'^'>p'\D^Pac){D^PM)} = 0, (2.18) 

(5) h[>fDlpab - e''p'%D^pac){D_pM) + e"{D_pab){D_a)] 

-^e'^^D+D^Pab + D^D+pab) + ^e''p'"^{iD^pac){D+pM) + iD^pbc)iD+Pad)} 
-^e'^{{D.p,b){D+a) + (D+p,,)(Z?_a)} 

+e'^(D-P„,)(Z?_/i) + ^e^^PacPbdF^^F^^ - \e^'' PabPcdF+lF^^ = 0. (2.19) 

Here R2 is the scalar curvature of A'^2) and the diffA'^2-covariant derivatives are given by|^, |lO|, 

F^l = d+A^ - d.A^ - [A+, AJ\l, (2.20) 

D^a = d±a-[A±, a]i^, (2.21) 

L»±/i = 9±/i- [A±, /i]l, (2.22) 

F)±Pab = d±Pab - [A±, p]Lab- (2.23) 

In general, the diffA'^2-covariant derivative of a tensor density fab— with weight w with respect to 
the diffeomorphisms of N2 is given by 

D±fab... = d±fab... - [A±, /]Lafe..., (2.24) 

where the bracket [A±, f]Lab- is the Lie derivative of fab- along A± := A^da, 

[A±, fUb- ■■= A^dJab- + fcb-daAl + fac-dbAl + ■■■ + w{d,Al)fab-. (2.25) 

For instance, the diffA'^2-covariant derivatives of the area element e"^ and the conformal metric pab 
which are scalar and tensor density with weight 1 and —1 with respect to the diffA'^2 transformations, 
respectively, are given by 

Die"" = d^e" - A^dce" - {daA^)(f, (2.26) 

F)±Pab = d±Pab - AldcPab " PcbdaAl - PacdbA^ + {dcA^)pab. (2.27) 

If one uses the Leibniz rule in (|2.26j) . then one has 

D±a = d±a - A^dcCT - daA^ 

= d±a-[A±, a]L, (2.28) 



which is just the equation ()2.21|) . 

The spacetime integral of the scalar curvature of the metric (\'2.1^ is given by 

/ = dudv cPy L + surface integrals, (2.29) 

where L is given by0, H, El E3, E3, Q 



■U'"^Pa,F^^F^'^+e'^{D+<j){D^<j) - 1 



e''' PabF+^F^l + e^D+a){D^a) - -e''p'''p'^{D+pac){D_pM) - e'' R2 



-2e^D^h){D^a) - he^D^af + he'^p'^'p^^D^pacXD^Pbd). (2.30) 

Each term in (|2.3fl)) is manifestly invariant under the diffeomorphisms of A''2, since the {y"'}- 
dependence of each term is completely hidden in the diffA''2-covariant derivatives. In this sense one 
may regard N2 as a kind of internal space as in Yang-Mills theory, with the infinite dimensional 
group of diffeomorphisms of A''2 as the associated gauge symmetry. Thus, the Einstein's gravi- 
tation of 3-|-l dimensional spacetimes is describable as 1-1-1 dimensional Yang- Mills type gauge 
theory interacting with 1-|-1 dimensional scalar fields a, h, and non-linear sigma fields pab whose 
interactions are dictated by the above Lagrangian density L. If one uses the diffA'^2 gauge freedom 
so that A^ = 0, then the metric 1)2. 1() becomes identical to the metric of the null hypersurface 
formalism studied in •^^. In this paper, however, I shall retain the A_!^ field, since its presence will 
make the coordinate choice less restrictive and the diffA^2-iiivariant Yang-Mills type gauge theory 
aspect more transparent. 

III. A SET OF QUASI-LOCAL BALANCE EQUATIONS 

Notice that the equations ()2.13j) . 1)2.14(1 and ((2. 17(1 are partial differential equations that are 
first-order in D_ derivatives. Therefore it is of particular interest to study these equations, since 
they are the analogs of the Einstein's constraint equations in the standard 3-1-1 formalism. Thus, in 
this (2,2) formalism, the natural vector field that defines the evolution is D_. Then the momentum 
variables vr/ = {7r/t,7r(j, vra,7r"^} conjugate to the configuration variables q^ = {h,a,A_^, pab} are 
defined as 

dL , , 



They are found to be 



T^h = -2e"(Ll-(j), (3.2) 

vr,, = -le^iD-h) - 2he''{D-a) + e^'iD+a), (3.3) 

vra = e2V,F+^, (3.4) 

^ab ^ he-pacpbd^j^p^^^ _ l^^pacpbd^jj^p^^y (35) 



Conversely, one can express D- derivatives of the configuration variables in terms of the conjugate 
momenta as follows, 

D-h = -^e-'vr, + ]^D+a + ^/le-'vr/,, (3.6) 

D^a = -]^Q-''^h, (3.7) 

F+_» = e-^V^^fe, (3.8) 

D-Pab = -f^e'" Pacpbd-J^""^ + ^D+Pab. (3.9) 

Notice that tt""^ is traceless 



Pab^"" = 0, (3.10) 



due to the identities 



p'^'D^pab = (3.11) 

which are direct consequences of the condition 

det pab = 1. (3.12) 

The Hamiltonian density Hq is given bvjlSl. Il4l| 

Ho := TTiD^q^ -L 

= i:f + total divergences, (3.13) 

where H is 

111 1 

+^7rj,{D^a) + l^T:-\D+pab) + ^e" p'^' p'\D+pa,){D+pM) + e''R2. (3.14) 

Notice that H and Hq are Hamiltonian densities that differ by total divergences only. In terms of 
the canonical variables, the first-order equations (|2.13jl . (|2.14|) . and (|2.17|) can be written as, after 
a little algebra. 



(i) TT^^D+pab + TTaD+a - hD+TTh " 8+ {h Tif, + 2e''D+a) 

+da{h-KhA^ + 2A^e"D+a + 2he-'' p^^^ + 2p'*^9b/i) = 0, (3.15) 

(ii) H - d+TTh + da{A_^TTh + e-'^p'^Kh) = 0, (3.16) 

(iii) d+TTa - db{A^TTa) - ndaA^ - TT^^aCr + daTTa - T^hdah - 'K^'^daPbc 

+56(7r^^p„,) + d,{TT''^Pab) - daiTT^'^Pbc) = 0. (3.17) 



Notice that (|3.15|) and (|3.16|1 are divergence- type equations of the following form|l5| 



A + d+B + daC'' = 0. (3.18) 
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One can also express ()3.17() as another divergence- type equation. If we contract H3.17() by an 
arbitrary function ^" of {v, y**} such that 

d+C = 0, (3.19) 

then we have 



7r'''£^Pab + 7Ta^a + 7Th£^h + 7T,£^A^-d+{e7ra)+da[-eTra + 27T''''epbc + A^eT^b) = 0, (3.20) 

which is in the same form as p. 18(1 . Here £^fab-- is the Lie derivative defined in ()2.25|) . 

£^fab- = [C, fhab- {^■=eda)- (3.21) 

Integrals of the divergence- type equations (|3.15|) . (|3.16j) . and (|3.20|) over a compact two-surface N2 
become, after normalizing by l/lGvr, 

d 



d_ 

du 
d_ 
du 
where we used the identities 



^f/(n, ^) = Y^ fd^y {TT^'D+pab + 7T^D+a - hD+iTh) , (3.22) 

P{u,v) = r^j(fyH, (3.23) 

L{U, V-i) = -^ jcPy (TT^^£^Pab + TTa£^<T " h£i:7rh - A_^£i:TTa) , (3.24) 



<i>(f'yKh£^h = — (p(fyh£^7rh, (3.25) 

j>(fyTia£(^A^ = -<j>d^yA^£^'Ka. (3.26) 
Here U{u,v), P{u,v), and L{u,v;^) are two-surface integrals defined as 

U{u, v):=-^j (fy (h Tih + 2(fD+a) + tj , (3.27) 

P{u,v):=^j>dFy{^h)+P, (3.28) 

Liu,v;0-=r^<fd^yie7ra)+L (9+^ = 0), (3.29) 
where U, P, and L are undetermined subtraction terms that satisfy the conditions 

du du du 
Notice that choices of subtraction terms are not unique, since it is the subtraction terms that 
define the references against which these quasi-local quantities are measured. A natural criterion 
for the "right" choice of subtraction terms would be that values of quasi-local quantities reproduce 
"standard" values in the well-known limiting cases, but otherwise they can be chosen arbitrarily. 
Let us notice that the integrand of the r.h.s. of 1)3.24(1 assumes the canonical form of angular 
momentum-fiux, 

Toar~E^^"^?9^ (3.31) 
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u = constant null 
hyp ersur face 




V = constant timelike 
hypersuxface 



V = constant spacelike 
hypersurface 



= constant null 
hyp ersur face 



(a) 2h>0 



(b) 2h < 



FIG. 2: (a) This figure shows the spacetime geometry in the region where 2h > 0, and f ==constant 
hypersurface is timehke. (b) This figure shows the spacetime geometry in the region where 2h < 0, and 
T;=constant hypersurface is spacehke. In both cases N2 is represented by a circle. 



where ^ is tangent to the two-surface A^2- One can also put the r.h.s. of ()3.22p into the canonical 
form of energy-flux, 






(3.32) 



where r/" = 5_^. To show this, let us contract H3.17|) with A^ and integrate over A^2- Then we 
have 



j(fy[A^d+-Ka) = <pd^yyn-''^£A+Pab + T^a£A+(^-h£A+TTh)- 



(3.33) 



If we use the equation (|3.33() and the diffA'^2-covariant derivative D^ defined in 1)2. 24() . then ()3.22|1 
can be written as 

Q^U{U, v) = Y^f d^y {-^""'d+pab + T^ad+a - hd+TTh - A^d+TTa) , (3.34) 



where the r.h.s. indeed assumes the canonical form of energy-fluxjl^. 

In the region where d+ is timelike (2/i > 0), the equations ()3.23() . H3.24() . and H3.34() are quasi- 
local balance equations that relate the instantaneous rates of changes of two-surface integrals at 
a given u-tirae to the associated net fluxes across the timelike tube generated by d+ (see FIG. 
2a). Let us remark that, unlike the Tamburino-Winicour's quasi-local conservation equations that 
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are "weak" conservation equations since the Ricci flat conditions (i.e. the fuh vacuum Einstein's 



equations) were assumed in their derivation |40[, our quasi-local balance equations are "strong" 



conservation equations in that only four Einstein's equations of the divergence-type were used in 
the derivation. 

In the region where 9+ is spacelike (2/i < 0), the vector field d/du is spacelike, 

< ^, ^ >= -2h + e'^pabA^A^ > 0, (3.35) 

so that u is the radial coordinate (see FIG. 2b). Then the equations ()3.15|) . ()3.16|) . and (|3.2U|) 
are constraint equations rather than balance equations, and each equation splits into a divergence 
term and a source term. The source term either assumes the canonical form of energy-momentum 
"density" T,TTi£^q and T,TTid^q as in Eqs. ()3.31|) and ()3.32|) . or is given by the Hamiltonian 
"density" H in 1)3.14(1 . Notice that the source term is not "flux" but "density", since the v = 
constant hypersurface is now a 3-dimensional spacelike hypersurface. These constraint equations 
describe how the quasi-local quantities at a given u-radius change as the radius u changes on a given 
spacelike hypersurface. Thus, the difference of two-surface integrals evaluated on two successive 
two-spheres on a given spacelike hypersurface is given by the 3-dimensional spatial integral of the 
"density" over the region between the two-spheres of interest. This is exactly what happens in 
Maxwell's theory, where the Gauss law constraint is given by 

V • ^ = 47rp. (3.36) 

If one integrates ()3.36|) over a 3-dimensional spacelike region V whose boundaries are two-spheres 
5i and ^2, then one has 

f En da — f En da = 4iT p dv. (3.37) 

JSi Js2 Jv 

Thus, in Maxwell's theory, the difference of two-surface integrals of the "momentum" En on two 
successive two-spheres is given by the spatial integral of the charge "density" over the volume 
between the two-spheres. In this paper, however, we shall be concerned with the case 2h > only, 
and the case 2/i < will be discussed elsewhere J4l|. 
If we introduce a function Wr{u,v) defined as 

Wr{u,v) := ^ J"du(j^d^yA^d+TTa, (3.38) 



then we have 



d_ 
du 



-Wr{u, v) = r^ jd^y A^d+^a, (3.39) 

so that ()3.34|) can be written as 

-[u{u,v)+WR{u,v)] = ^jd^y{Ti'''d+pab + TTad+a-hd+7ih). (3.40) 



d_ 
du 
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-fwniu, v)r^Y. ^"d+La, (3.41) 



Notice that the r.h.s. of ()3.39|) has the form 

d_ 
du 

where 

9." ~ A^, (3.42) 

La ~ TT^vr^, (3.43) 

which represents the work done per unit n-time by changing the angular momentum La of the 
system that has the angular velocity fi". From this perspective, the equation (|3.39j) is the work 
done on A'^2 per unit u-time by changing the angular momentum density VTa/lSvr of gravitational 
field that has the angular velocity A^ at each point of A'^2. This observation suggests that Wr{u, v) 
be identified as the quasi-local rotational energy of A^2. Indeed, as is shown in section IVll Wfi{u, v) 
reduces to the Carter's constant for the asymptotically Kerr spacetimes, the total angular momen- 
tum squared [23, |2a, 1221, l2a|. This is a strong indication that supports our identification of Wr{u, v) 



as the rotational energy of gravitational field in the circumstances where no isometrics are present. 
In the limit where A^ is independent of u such that 

d+A_^ = 0, (3.44) 

Wr{u,v) becomes 

WRiu,v) = -^<f d^y (A^TTa) - -^ (f Sy{A^TTa). (3.45) 

iDVr Ju=u lOTT Ju=Q 

IV. HAMILTON'S EQUATIONS OF MOTION 

Let us define the Hamiltonian K as the integral of H in 1)3. 14(1 . 

K:= jdujd^y[H + \{detpab-l)}, (4.1) 

where A is a Lagrange multiplier that enforces the unimodular condition (|3.12l) . In the Appendix, 
we have shown that the equations (|2.15() . (|2.16|) . (|2.18l) . and H2.19() are the Hamilton's equations 
of motion 

-'' - £• '-) 

5K 
D_„ = -^, (4.3) 

where {vrj, q^} are 

vr/ = {tt/,, tt^, TTa, tt"^}, q^ = {/i, a, A^.pah}, (4.4) 
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assuming the boundary conditions 

Sa = Spab = (4.5) 

at the endpoints of the n-integration. Thus, together with the divergence-type Einstein's equations 
(tm . (tmi) . and (tTTTl) . from which follow the integral equations (TT^ . (IT^ . and (TT^ that 
may be interpreted as either the quasi-local balance equations or constraint equations depending 
on the signature of the 3-dimensional hypersuface, the Hamilton's equations of motion H4.2|) and 
H4.3|) make up for the complete set of the vacuum Einstein's equations. Thus, in the (2,2) fibre 
bundle formalism, the Einstein's equations split into twelve first-order Hamilton's equations of 
motion dictating the evolution along the out-going null geodesic and the four quasi-local balance 
equations or the constraint equations that implement the Hamilton's evolution equations. 

V. GEOMETRICAL INTERPRETATIONS 

Two-surface integrals ()3.27() . (|3.28() . and H3.29() can be given geometric interpretations in terms 
of the area of A''2 and null vector fields orthogonal to A''2. In order to see this, it is necessary to 
recall the definitions of the in- and out-going null vector fields n and / given by (|2.7j) and (|2.8j) . 
respectively. 

A. Quasi-local energy 

Let us observe that, apart from the subtraction term U, 1)3. 27(1 can be written as the Lie 
derivative of the area A of A''2 along n. Notice that we have 

i(fy(hirh + 2e''D+a^ = 2 i cfye" [o+a - hD_a) 

= 2<fd^y£ne\ (5.1) 

But one has 

j cPy £n(f = £nA, (5.2) 

where A is given by 

A= ixfye". (5.3) 

The identity 1)5. 2|) follows from the fact that the order of tf'y integration and the Lie derivation £„ 
is interchangeable, since the in-going null vector field n is orthogonal to A'^2- Thus we have 

U{u,v) = ^£nA + U. (5.4) 

ovr 

In order to fix [/, it is necessary to introduce a reference spacetime. In principle, the reference 
spacetime can be chosen arbitrarily, provided that the pull-back of the background metric to A'^2 
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is the same as e" pab- If we denote the coordinates of the reference spacetime as (u, v, y"), then its 
metric can be written as 

df' = -2dudv - 2hdu^ + e^Pab {dy" + A_^du + A_^dv) (dy^ + A^du + A}^dv^ , (5.5) 

where {h,A_^} are the embedding degrees of freedom of A'^2 into the reference spacetime. The 
vector fields {n, 1} 

are null with respect to the background metric, and satisfy the same normalization conditions as 
before, 

< n, n >ref= 0, < /", r>ref=0, < n, r>ref=-l. (5.8) 

If U is chosen as 

U := -^£nA (5.9) 

Sir 

such that it satisfies the ^-independent condition 

dtJ 

- = 0. (510) 



then (|5.4j) becomes 



U{U,V) = —£(_ri-n)-^- (5-11) 



■iTT 



This expression is entirely geometrical, stating that U{u, v) is determined by the rate of change of 
the area A of N2 along the difference n — n of the in-going null geodesies, and becomes zero when 

n = n. (5.12) 



B. Quasi- local linear momentum 

One can also express P{u, v) geometrically. Let us notice that 

Therefore, if we choose the subtraction term P as 

P := —£tA (5.14) 
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such that it satisfies the u-independent condition 



AP = 0, (5.15) 



then ()3.28|1 becomes 



P{u,v) = -^£f^i_j)A. (5.16) 

Thus, P{u, v) is given by the minus of the rate of change of the area A of N2 along the difference 
/ — / of the out-going null geodesies, and becomes zero when 

I = I. (5.17) 



C. Quasi-local angular momentum 

Let us now find the geometrical expression of L{u,v]£,). If we notice that the Lie bracket of n 
and / is given by 

[n, /]l = -F^^da + {D^h)l, (5.18) 

then we have 

= --^fd'ye^Un, Tl 

where we used 

Ca = e'^Pab^', (5.20) 

and the fact that / is orthogonal to S,, 

na = 0. (5.21) 

If we choose the subtraction term L as 

^■=^ jd^y^" < t [n, 1]l >ref, (5.22) 



and require that the condition 



^1 = (5.23) 



hold, then (|3.29j) becomes 

L{u,v;C) = —^<f(fye''(^<C, [n, 1]l > - < t [n, I]L>rc{). (5.24) 
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Thus, L(n, v; ^) is given by the integral over A'^2 of the projection of the twist [n, /]l onto ^ = C"'da 
modulo the subtraction term L. Notice that ^" is an arbitrary function of {u,y^}, so that the 
vector field ^ need not satisfy any Killing's equations. Thus, ^ is an arbitrary vector field tangent 
to N2, defining the direction of rotation at each point of that surface. 

If L{u,v;S,) is to be regarded as an acceptable candidate of quasi-local angular momentum, its 
value must be zero for any two-surface embedded in the flat Minkowski spacetime. Our expression 
clearly satisfies this criterion, as can be seen by the following observation 42]. Let be a diffeo- 
morphism from a given spacetime M3+1 to itself, and 0* be the push-forward associated with (j). 
For any vector fields X, Y defined on Ms+i, the inner product and the Lie bracket are preserved 
by the mapping (f)^ 4j|, 

< (p^X, (P^Y >= (/>, < X, y >, (5.25) 

[<P,X, <P,Y]i^ = <P4X,Y]i^. (5.26) 

Suppose that M3+1 is the flat Minkowski spacetime, and let n and I be the null vector fields at 
each point of A'^2 of the flat Minkowski spacetime. Then these null vector fields remain null and 
the normalization condition is preserved under the mapping 0*, since we have 

< (/)^,n, (/)*n >= (j)^ < n, n >= 0, (5-27) 

< 4>J, (pj >= (p^ <l, I >= 0, (5.28) 

< 4>^n, (pj >= (p^ <n, I >= -1. (5.29) 



The following identity 



6*n, (pj]i, = (p^,[n, /]l = (5.30) 



is also true, since two null vector fields at any points in the flat Minkowski spacetime must commute. 
Therefore, [n, /]l = holds on any two-surface (and its deformations) in the flat Minkowski 
spacetime. Thus, L{u,v;S,) is zero on any two-surface N2 in the flat Minkowski spacetime, modulo 
the subtraction term that can be trivially put to zero. 

Notice that L(^) is linear in ^. That is, for any S, given by 

e = <! + 6^2, (5.31) 

where a, b are constants, we have 

L(e) = aL(6)+6L(e2), (5.32) 

which shows that the quasi-local angular momentum is additive. Asymptotic properties of quasi- 
local angular momentum and its flux will be studied in section IVll 
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D. Quasi-local rotational energy 

Like quasi-local angular momentum, the value of any reasonable candidate of quasi-local rota- 
tional energy must be zero for any two-surface in the flat Minkowski spacetime. The quasi-local 
rotational energy Wr defined in (|3..S8j) trivially satisfies this criterion. Since Wr can be written as 

Wr = -r^ fduj<fy A^d+if'' Pab[n, l]i), (5.33) 

it is zero when 

[n, l]l = 0, (5.34) 

which is true for any n and / of the flat Minkowski spacetime. 

VI. ASYMPTOTICALLY FLAT LIMITS 

In this section I will show that the limits of quasi-local balance equations in the asymptot- 
ically flat zones are the Bondi energy-loss, linear momentum-loss, and angular momentum-loss 
eauations|13l. Il4| . Moreover, the integral Wr turns out to be proportional to the total angular mo- 
mentum squared in this limit, which is a strong indication that it is to be regarded as a quasi-local 
generalization of the Carter's constant |2fil. Yin. l27| for a generic gravitational field. 



The general form of asymptotically flat metrics|l3l. Il9l. l3ll. l32l. l33l| is given by 



:or a ge 



, r, , , / 2?7i \ r, /4r?T,asin^t? Ama? sm^'dcos'^'d \ , , 

ds^ — > -2dudv - 1 \ ]du^ + { \ ]dudip 



^y2[l + l^^ + ...y^2^^2^:^^2^^^^<^^,__y^2 



+sin 



9 -I )"T 

y^ J \ y 

'Ama^ Sm?a^ \ , , /a^sin^i? 



„/^rna am a \ , , /a sm w \ , o , . 

K-^^ + -;;^ + ---yvd^-[-^^ + ---)dv^ (6.1) 

as t; ^ oo. From this expansion, asymptotic fall-off rates of the metric coefficients are found to be 



e"" =t;2sim?|l + 0(^)}, (6.2) 

p,, = (sin^)-i{l + £Ml^ + o(J_)\ (6.3) 

L y y^ J 

P,, = sin^{l-^(^^lA^ + 0(J,)}, (6.4) 

L y y^ J 



V 

P% = C'(^), (6.5) 

Im 1 

2/1 = 1 - + 0{^), (6.6) 

^/ = ^ + 0(^), (6.7) 

2ma^ ^, 1 



^_^^ ^^ (6.8) 

^/ = 0(4), (6.9) 



•^6' 
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and their derivatives are given by 



a+a = o(-), 

9 1 

d.a = - + 0{^ 



d+Pa 



o(i), 

1 



yZ 

£iPab = O(-), 



vr 



1 



^e>V^(9+Pcd) + 0(l), 
-4T;sim? + 0(l), 
-2t;sinT? + 0(l), 



7r<o = Gmasin^?? + 0{ 



1 



1, 



^^ = 0(^). 



Therefore, n and / become, asymptotically, 



d 



1 m\ 9 



n 



a /i m\ a 
(9?y, V 2 ?; / (9?; ' 



9^1 
The natural reference spacetime at the asymptotic infinity is the flat Minkowski spacetime, 



dr 



-2dudv - du^ + r^{d-d'^ + srn^MLp^), 



(6.10) 

(6.11) 

(6.12) 

(6.13) 

(6.14) 

(6.15) 
(6.16) 
(6.17) 
(6.18) 

(6.19) 

(6.20) 
(6.21) 

(6.22) 



where N2 = 82- Thus the embedding degrees of freedom of 5*2 into the flat Minkowski spacetime 
are given by 



so that n and I are given by 



:±" = 0, 


2h = 


d 


1 d 


n = — 


— 


du 


2dv 


dv' 







(6.23) 

(6.24) 
(6.25) 



A. The Bondi energy-loss relation 

In the asymptotic region ?; — > 00 where v = constant hypersurface is timelike, the r.h.s. of 
1)3.34(1 represents the canonical energy- flux carried by gravitational radiation crossing A'^2- Then 
the l.h.s. should be identified as the instantaneous rate of change in the gravitational energy of 
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the region enclosed by N2. Energy- flux in general does not have a definite sign, since it includes 
the energy- flux carried by the in-coming as well as the out-going radiation across A'^2- In the 
asymptotically flat region, however, energy-flux turns out to be negative-definite, representing the 
physical situation that there is no in-coming flux coming from the infinity |43|- 

Let us show that the balance equation H3.34() indeed reduces to the Bondi energy-loss formula 
at the null infinity. Let U be given by (|5.9)) and use n in (|(i.24j) . Then we find that 

U = l. (6.26) 

Let us suppose that the coordinates {u, v, y"} approach the coordinates of the Minkowski spacetime 

{n, u, ■(?, (/?} as f — > 00, 

u — >u, V — >v, y" — > {^, ip]. (6.27) 

Then U trivially satisfies the condition (|5.10j) . since 

^-U ^^U = Q. (6.28) 



If we use the asymptotic formula 



then the total energy is given by 



du du 



m d , ^ 

^-^^-TT.^ (6-29) 



V 



dv'' 



Ub{u) := lim U{u,v) = lim -—£,n-n)-^ 

= m(n), (6.30) 

which is just the Bondi mass at the null infinity. 

Asymptotic limit of the balance equation 1)3.34(1 is found to be 

-^[/B(n) = - hm -i- /d(7^VV^(9+P6c)(5+Pad) 

= - hin :^idn v\j^%U\) < 0, (6.31) 

where j_^jj is the shear current defined as 

U\:=p-'^d+p,, (iA = 0), (6.32) 

which represents traceless shear degrees of freedom of gravitational radiation. The relation (|6.31|) 
is just the Bondi energy- loss formula with the correct normalization coefficient |43]. Notice that the 
energy-flux is bilinear in j_j"j. Equivalently, it can be written as 

^UB{u) = --^idn{d+Cf<0, (6.33) 

du Ibvr Js2 

if one uses the expansion of pab given by (|6.3|) , (|6.4|) , and (|6.5j) . 
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B. The Bondi linear momentum and linear momentum-flux 

The total linear momentum Pb{u) and its flux are trivially zero, 

^£ - 



PbH := lim P(n, t>) = - lim —fnnA = 0, (6.34) 



^PB(n) = 0, (6.35) 

since we have 

I -I — >0. (6.36) 

That the net-flux of the total linear momentum is zero can be also seen by evaluating each term 
of H in ()3.14l) in the asymptotic limit. Let us notice that although the fourth, sixth, and seventh 
term in (|3.14l) are not zero individually, they add up to zero asymptotically, 

^e-V,,p,,vr-^^^ + l-^-^iD+p^,) + i-e>"V^"(^+Pac)(I?+P6d) 
= ^eV'^V^'(^-Pac)(I)-P6d) 

= O(^), (6.37) 

where we used the definition of vr"'' in 1)3. 5(1 . The third and fifth terms become zero, respectively, 

_le-V'vra^6 = 0(^), (6.38) 

^7Th{D+a) = 0{^). (6.39) 

The remaining non-vanishing terms are given by 



lim --^ / d^y (le-'^vr^vr,) = -1, (6.40) 

1 /* 1 1 

lim i Sy f-Zie"'"^?') = -, (6.41) 

lim -^<fcP ye'' R2 = lx, (6.42) 

where x = 2 for S'2. Since these terms add up to zero, it follows that the net momentum-flux H is 
zero at the null infinity. 

C. The Bondi angular momentum and angular momentum-flux 

Likewise, the total angular momentum L-q{u] ^) is defined as the asymptotic limit of quasi- local 
angular momentum L{u,v;S,), 

Lb{u;0 '■= lim L{u,v;i) 

1 ^2_. 



■ Jkn — ^d"ye"( < i, [n, /]l > - < i, [n, /]l >rcf )• (6.43) 
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Let ^ be asymptotic to the aziniuthal Killing vector field 



( = ea.^-^. (6.44) 



From the expansions (|6.2|) . • • •, (|(i.9|) . we find that 



e'' — > u^sim?, (6.45) 

C^ — > w^sin^T?, (6.46) 

C^ ^ 0, (6.47) 



where we used the definition of ^a in (|5.2()j) . The Lie bracket [n, /]l is found to be 

'"■'K--^ + o(;^). (^■^«) 



[n, /]^^0(1). (6.49) 



Since [n, /]l = 0, we find that 

c27r rn 



1 rZTT PTT 

Lb{u;C) = dip di? (6ma) sin^i? 

IQlT Jo Jo 



10 JO 

= ma, (6.50) 

which is just the total angular momentum of the Kerr spacetime. 

The Bondi angular momentum-loss relation will be obtained by taking the limit of the quasi- 
local balance equation 1)3. 24() . 

^ = Jn^ -^ j (fy {Tr'''£ ^pab + T^^^a - h£^TTh - A_^£^7ra) . (6.51) 

Let us evaluate each term of this equation. The first term has a finite limit, which is 

icPyTT'^'X^Pab ^-\id^ v^p'''p'\d+p,d){£^Pab), (6.52) 

JS2 ^ JS2 

where we used the notation 

£<fi '■= £d/dip- (6.53) 

Notice that 

TT^i^^o- — > {-2vsm^ + 0{l)}£^a. (6.54) 

But a becomes, asymptotically, 

a — >2lnv + ln\sm'd\ + \nh+0(^'^}, (6.55) 

so that we have 

^v=^ = 0(^)- (6.56) 
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Thus the second term becomes zero, 



i(fy-K^£^a = 0(^-^ — ^ 0. 



Notice also that 



h£^TTh 



£^{h7Th) -■Kh£^h 
— > £^{h7rh) — 4:sm'd£^m + 0(-] 

= £^{h7rh-^msm-d) +o(-^ 

Thus the third term also becomes zero, 

j<fyh£^^h = 0{^) ^0. 

The fourth term dies off much faster, since 



J S2 ^ 



0. 



From (llOa . (ROTTl) . dniMI), and (IfO)!]!) . we find that (jH3T]) becomes 



dLv 






(6.57) 



(6.58) 



(6.59) 



(6.60) 



(6.61) 



which is just the Bondi angular momentum-loss relation with the correct normalization coefficient. 
It is worth notin g th at 1)6. 61() is the coordinate-dependent expression of the angular momentum-flux 
discussed in 
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231 . |24 1 . If we use the asymptotic expansion of pab given by (|6.3(1 , 1)6. 4() , and (|6.5() , 
then this relation can be expressed as 

dLB _ 
du iDvr J S2 



--^ i dn {d+c){£^c). 

167r J So 



(6.62) 



D. Gravitational Carter's constant 



Let us find what Wr becomes in this limit. The equation (|3.39|) becomes 

3 



du 



„ , dip di} siir^ [ma) —^ — - 
Attv^ Jo Jo du 

'^ d 2 

v-^ du 



If we choose the constant of -u-integration as zero, then we have 



lim u VFr 

D— »00 



[ma] 



(6.63) 



(6.64) 



which is just the total angular momentum squared for the Kerr spacetime. Thus, Wr may be 
regarded as a quasi-local generalization of the Carter's constant 25|, |2g, |27|, |2^, and physically, 
could be interpreted as gravitational contribution to the quasi-local rotational energy. 
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VII. RELATION BETWEEN ENERGY-LOSS AND ANGULAR MOMENTUM-LOSS 

In general, if a given system undergoes an energy-losing process, then it always accompa- 
nies angular momentum-loss, unless the system remains spherically symmetric throughout the 
whole process. Since we already found general expressions of quasi-local energy-flux and angular 
momentum-flux, it is natural to ask what relation exists between them, if there is any. The relation 
is that the angular momentum-flux 1)3. 24() and energy-flux (|3.34() transform into each other 

^L{u,v;O^^U{u,v), (7.1) 



under the exchange of the derivatives 



^« " I; P-^' 



in the flux integrals. 

Let us examine the implication of this symmetry for spacetimes close to a background spacetime 
that possesses two commuting Killing vector fields. Let us choose the coordinates of the background 
spacetime as {u,v,'d,if}, and let {d/du,d/dip} be two Killing vector fields of the background 
spacetime. Let us also suppose that the coordinates {u,v,y°'} approach the coordinates of the 
background spacetime {u,v,'&,ip} as v — > oo, 

^^^2, v^v, y''^{'&,^}. (7.3) 

If we perturb this background spacetime by adding a small amount of gravitational waves, then 
we may regard these waves as propagating in the background spacetime, carrying a small amount 
of energy and angular momentum. Let us write q^ = {h,a,A_l^,pab} and ttj = {vr/^, tTq-, vr(j,7r°''} 
about an exact solution {q^ ,Tti} of the Einstein's equations, 

q' = q'{v,^) + 6q', (7.4) 

TTj = 7ri{v,'d) + 67ri, (7.5) 

where {dq^,67:j} represents gravitational waves propagating on the background spacetime. The 
dependence of a given mode on u and ip is given by 

5q^ = Q\v, z?) 6^'""+™^^ + c.c, (7.6) 

6^1 = n/(v, 1?) 6*'^"+*™^^ + c.c, (7.7) 

where {Q , 11/ } are the amplitudes of the wave that has the frequency uj and the azimuthal quantum 
number rriz [rriz = 0, itl,---). Now, if we use (|7.3jl and the fact that {d/du,d/dip} are the 
timelike and azimuthal Killing vector fields of the background spacetime, respectively, then from 
the symmetry 1)7. If) and 1)7. 2() we obtain the following relation |4J] 

dU UJ dLz 
du THz du 
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which is a wen-known relation between energy-loss and angular momentum-loss for perturbations 
around the stationary and axi-symmetric spacetime. Thus, the symmetry ()7.1() and (|7.2() is the 
sought-for relation between the energy-loss and angular momentum-loss for a generic gravitational 
radiation that has no isometrics. 



VIII. DISCUSSIONS 

The key ingredient of the (2,2) fibre bundle formalism discussed so far is the observation that 
the out-going null vector field defines a natural time flow. With the affine parameter of the out- 
going null geodesic as the time function of the theory, the canonical variables were introduced 
and the non-vanishing gravitational Hamiltonian was spelled out. Then I obtained the Hamilton's 
equations of motion from the Hamiltonian by the variational principle, which are the evolution 
equations along the out-going null geodesies with respect to the affine parameter. Thus, in this 
paper, I have shown that the Einstein's equations split into twelve first-order Hamilton's equations 
of motion and the four quasi-local balance equations or constraint equations that implement the 
Hamilton's evolution equations. 

I also found coordinate-independent and geometric expressions of quasi-local gravitational en- 
ergy, linear momentum, and angular momentum for any two-surface. The corresponding fluxes of 
gravitational field were found to assume the canonical form of energy-momentum flux, 

rOa??"~5]^/£y, (8.1) 

/ 

just as in standard field theories. I have shown that the quasi- local balance equations correctly 
reproduce the well-known Bondi relations at the null infinity of asymptotically flat spacetimes. 
However, because of the breakdown of the coordinate system due to potential occurrence of caustics 
after a finite propagation along the out-going null geodesies, there could be difficulties in extending 
the Hamilton's evolution equations beyond the caustics. In principle, however, it is still possible 
to approach the null infinity by using a new coordinate system after the breakdown of the old 
one, and perhaps one could use the quasi-local balance equations across the caustics and search 
for "weak" solutions 45|. However, it must be mentioned that, the farther out one goes, the less 
likely is the chance for caustics to occur due to the weakness of gravity near the infinity. If one 
is interested in the strong gravity region near black holes, or black hole dynamics itself, then the 
caustics might cause serious problems since they are much more likely to occur as we approach 
strong gravity region along the in- going null geodesies. 

Quasi-local angular momentum was defined in this paper for spacetimes that have no isometrics, 
and was found to be zero for any two-surface in the flat Minkowski spacetime. It was found to 
have the additive property, being a linear functional of a vector field ^ that defines the rotation 
at each point of the two-surface. One might be interested in studying symmetry properties of this 
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quasi-local angular momentum, and look for some generalization of the BMS symmetries at a finite 
region [l9|. 

In addition, I obtained a quasi-local generalization of the Carter's constant of gravitational field, 
and interpreted it as gravitational contribution to the quasi-local rotational energy. The Carter's 
constant is known to exist when the system under study has two commuting Killing vector fields, 
such as the Einstein-Maxwell system and the Einstein's equations coupled to a scalar field. For 
a generic gravitational field that has no isometrics, no analog of the Carter's constant is known. 
In this paper, I presented a candidate for the generalized Carter's constant which becomes zero 
for any two-surface in the flat Minkowski spacetime, and reduces to the total angular momentum 
squared in the asymptotic region of the Kerr spacetime. It is interesting to see how this quasi- local 
Carter's constant generalizes when the electromagnetic and scalar fields are present. 

There could be a number of applications of the quasi-local balance equations in astrophysics. 
The most important and challenging problem seems to be the calculation of back-reaction on the 
geometry of black holes as a consequence of the emission or absorption of gravitational radiation. 
One could also use the quasi-local balance equations in searching for consistent boundary conditions 
at a finite boundary in numerical relativity, since the boundary data at a finite boundary must 
satisfy the quasi-local balance equations. These problems are left for future works. 

Another issue in this (2,2) formalism is the well-posedness of the initial value problem. When 
the initial 3 dimensional hypersurface is chosen spacelike, there is no problem in the well-posedness 
of the initial value problem since the null direction can be viewed as a limit of timelike direction. 
But there are several other choices of initial surfaces, such as the double null initial surfaces and 
the initial/boundary value problems where the boundary can be either timelike or null. One of the 
difficulties associated with the characteristic or initial/boundary value problem is that one has to 
know the "flows of information" across the characteristic or the timelike boundary that belongs 
to the future. In these hybrid formulations of the Einstein's equations, not so many articles that 
aim to study the well-posedness of the field equations appeared. However, a series of the recent 
papers by Frittelli|46|| shows that, for a certain choice of first-order variables for the characteristic 
problem of the linearized Einstein's equations, the system can be cast into manifestly well-posed 
form. For the non-linear characteristic problems, the notion of well-posedness is still not available. 
It is interesting to examine whether the first-order variables in this paper might have any relevance 
in establishing the well-posedness of the non-linear characteristic initial value problem. 

Finally, there are problems related to the gauge invariance of this (2,2) fibre bundle formalism. 
It is obvious that this formalism is tied to a particular gauge, and the non-vanishing Hamiltonian 
is obtained as a consequence of selecting a particular time function, namely, choosing the affine 
parameter along the null direction as the time function. But one should notice that, in the standard 
ADM formalism, it is also possible to obtain another non-vanishing Hamiltonian if one chooses a 



time function such as the Gauss normal time coordinate 29|. Moreover, if one quantizes the theory 
in a particular gauge, the resulting quantum theory will depend on that gauge, losing the spacetime 
diffeomorphism invariance that one wishes to carry over to the quantum regime. In view of the 
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present situation that there is not any single complete version of sensible quantum theory of gravity, 
however, this gauge problem does not seem to be an urgent problem. Clearly, quantizing the full 
Einstein's gravity is beyond the scope of the present paper. 

Ackno^vledgments 

It is a great pleasure to thank S. Hollands, R. Geroch for interesting discussions at the last 
stage of preparing this paper, and in particular, R. Wald for the hospitality extended to the author 
during the sabbatical leave at the University of Chicago. The author also thanks the referee for 
his criticism that helped to improve the original manuscript on several delicate issues. 



28 



APPENDIX: HAMILTON'S EQUATIONS OF MOTION 

In this Appendix, I will show that, the Einstein's equations, ((TT3|) . ((TT^ . (|nH|l . and (fTTH)) . 
which are second-order in D_ derivatives, are the Hamilton's equations of motion, 

---^^ 

if the boundary conditions 

5a = 5pab = (3) 

are assumed at the endpoints of u-integration in K, where the Hamiltonian K is given by 

K = jdu jdPy [h + \ (det pab - 1) } , (4) 

and the Hamiltonian density H is 

+i7r,(D+cT) + ^T:-\D+pab) + ^eV'^V'^'(^+Pac)(I)+P6d) + e'^ii2. (5) 

(i) Variations with respect to vr/i and h 

It is trivial to see that the equation 



is identical to the equation H3.6|) . and the equation 



-- - -^ (^) 



can be written as 



Let us show that the equation Q is just the equation H2.15() . using the equation ©. Notice that 
each term in H2.15() becomes 

(i) 2e''D^_a = -^e-''7:l-D.7rh, (9) 

(ii) e^D^af = \e-'^7rl (10) 

(iii) ^e''p-'p'''{D.pac){D^pbd) = ^e-'^PabPcdT^'^'Tr"^ + ^TT^'D^p^b 

+^ey-'p'''{D+pac){D+pb,). (11) 
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From @, (fTU)) . and (fTT|) . it follows that the equation H2.15() is identical to the equation ©. 
(ii) Variations with respect to tTo- and a 
It is trivial to show that the equation 

is identical to the equation (|3.7p . In the variation 

I?-^. = -^, (13) 

the less trivial part is the following one, 

5 duwd y{-KhDj^a) = duid y7rhD+6a 

= — duxd'^y {D^TTh)5a + duxd^y Dj^{TTfi5a) 
= - du(td^y{D+Trh)5a + jdu—i id^yvr/i^o-j. (14) 

Therefore, if we assume the boundary condition 

<5a = (15) 

at the endpoints of u-integration, then we have 

^-^[ duj>d^yiThD+cj^ = --D+7rh. (16) 

The remaining variations are straightforward, so that ()13() becomes 

+iD+7r, - ^e"p-'p'\D+pac){D+p,a). (17) 

In order to show that the equation H17() is the same as the equation ()2.18|) . we need to express the 
derivatives of D- and D^ in H2.18() . using the conjugate momenta. Notice that the first term in 
((TTH|) becomes 

= -e^'^TTh-JTa- D^-K^ + e''D^D+a + -TThD+a + he~''7rl + hD_Trh- (18) 

Since the third term in the r.h.s. of ()18() can be written as 

= -^D+TTh + ^TThD+a + daie-y^'TTb), (19) 
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(|TH|) becomes 

(i) 2e"Dlh = -D^TT^ - e-'^TThTT, - -D+tth + nnD+a + -/le-'vr^ + ^e->afe/Ocd^"'vr'"^ 

+l.Tr'^''D+pab + |^e>'^V''(^+Pac)(Z?+PM) + a,(e->'^^6), (20) 

where we used the equation of motion of TTh given by 0- It is straightforward to express the 
remaining terms in (|2.18|) using the canonical variables. They are given by 

(ii) 2e^D^h){D^a) = ^e~^7Th7T^ - \^hD+a - ^he-'^nl (21) 

(iii) e^D+D^a = -^D+TTh + ^TThD+a, (22) 

(iv) e^D^D+a = -^D+TTf, + ^TT^D+a + da{e-^ p^'n), (23) 

(v) e^D+a){D-a) = -^TThD+a, (24) 

(vi) ieV'^V^'(^+Pac)(I?-P6d) = ^^"'D+pat + ^eV«V'^'(^+Pac)(I?+P6d), (25) 

(vii) e^'' PabF+IF^'L = e^2-pa^^^,. (26) 

If we plug (EOl), • • •, (Eni) into ((TTH)) . then the equation (|nH)) becomes 

+^e>V'(^+Pac)(I?+P6d) - 2/ie-{(i?2 ^) + i(Z)_a)2 + ipV'(i^-Pac)(i?-P6d)} 

= 0. (27) 

But the term in the bracket {} in l\27\i is zero if we use ()2.15p . and this shows that the equations 
((TTH)) and (HH) are identical. 

(iii) Variations with respect to iTa and A^ 

The equation 

5K 
D^A^ = ^ (28) 

on a 

is 

D^A^ = -e-^"" p^^TTb, (29) 

which is the defining equation (|3.8|) of tTq, since the covariant derivative of A_^ is given by 

D_A^ := F_^. (30) 

In order to write down the equation 

«-.. = -^, (31) 
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one needs to do the following variations. Notice that 

5Jdu<j^(fy (TThD+a) = Jduj^d^y {iThS^d+a - A^daO - daA^) } 

= jdujd^y { ( - -KhdaCT + da^h) 5A^ - da (j^hSA^) } . 

Thus we have 

A /* /* 1 1 1 

— — duid'^y-{iThD+a) = --nhdaCr + -daiTh- 

Let us also notice that, for an arbitrary field <^"'', the following is true, 

jdujd^C'^KD+Pab) 

= jdujd'^y C''[ - {5A^)dcPab - {dJA^)Pcb - {di,5A^)pac + {dc5A^)pab] 

= jdujd^y [{ - C'd.pab + daiC'pcb) + d,{C'pac) - 5e(CVa6)}M+^ 

-da[C'pcb^A^)-d,[c'pac5A^)+d,[c'pab5A, 

Let us consider the following two cases, (a) and (b). 
(a) If we choose 



(32) 



(33) 



(34) 



/-ab , „ab 

^ -=2^" ' 



then from 1)34^ we obtain 
6 



5A^ 



- Jdufd^y ^-K^^D+p,,) = -l-n'^daPbc + db^TT^'^f 



where we used the identity. 



.vr"^ = 0. 



(35) 



(36) 



(37) 



(b) If we choose 



^ab .^ l.^^paepbf^j^^p^^^^ 



which now depends on A_^, then 1)34^ becomes 

5 



SA 

] 
~4h 



- jdujd^y ^e'^p''^p'%D+p,,){D+p 
i-eV^'^p-(Z?+p,e)(5aP6c) + afe(i-eV^=D 



+ Pcaj, 



where we used the identity 



(38) 



(39) 



p'^'^D+Pab = 0. 



(40) 
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From the equations (|5H|) . (|H^. and (fS^ . we find that the equation (|5T|) becomes 

-db{\K''Pca + ^e>^^(l)+pea)}. (41) 

Using the definitions of the momenta ()3.6|) , • • • , (|3.9|) , one can easily show that the equation ()2.16|1 
is identical to (|41() . 

(iv) Variations with respect to vr and p^fe 

It is trivial to see that the equation 

^-P«^ = ^ (42) 

is just the equation ()3.9() . Let us show that the equation 

D.tt"'' = — ^ (43) 

is identical to the equation H2.19() . If we vary terms in K which do not contain D^pab, then we 
have 

du(bdy\^--e p TTciTd + ^e PcePdfT^ 7^2 7rc7rd + -e Pcd^r tt . (44) 

Varying terms linear in D^pat, we find that 

6jduj^d^y (i-TT'^^I^+p^,) = - lduj^d^yD+[^7r'''')5pab + jdu^[jd^y (i-vr'^^^p.fe) }. (45) 
If we assume the boundary condition 

5pab = 0, (46) 

then we have 

^/..„/d^.(i;,-D,,..)=-0,(±.-). (47) 

Now let us define 

S\ := p'^'D+p.b. (48) 

Then we have 

i-e>'^V^'(^+Pac)(I?+P6d) = ^e'^S^a S\, (49) 

and the variation of S\ is given by 

5S\ = -p'''S%5p,d + p'^'D+Sp.b. (50) 
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Therefore, we have 



jdujify { - ^^"p^'Si S\ 5pab + ^e''p''''S\D^5pab} 

jdujd\{- ^^""P^'Si S\-D+[^e'^p-^S\)]5pa,. (51) 



Therefore, we find that 
5 



6p, 



ab 



^du<fd^y{jj-e'^p'^'p''fiD+p,,)iD+p,f)} = -j^e'^ p^'^p"'p'fiD+p,,){D+paf) 



-Dj^e- 



V4^e>-p'"^D+pedj. (52) 

Finally, we have to vary the Lagrange multiplier term in (|4.H) . It is given by 

^ Jdu^fd^y A (det p,rf - 1) = A Z''. (53) 

The scalar curvature term e'^i?2 is a topological density that does not contribute to the metric 
variation. Prom (jlH), (gZI, ^^, and ((SHI), we have 

1 1 11 

D^TT-' + -e-^V-p^^vTevr, + -e->,,^-7r'"^ - ^+{^vr'^* + —e^p-^p^^D+p^)} 

-ll^y'"'p'"'p'HD+Pce){D+pdf) + \p'' = 0. (54) 

The Lagrange multiplier A is determined by taking the trace of (|54|) . Notice that for any traceless 
field x°'^ such that 

PabX"" = 0, (55) 



we have 



Thus, for x"''' defined as 



PabD^x"" = -X^'D^pab. (56) 



one has 



X"' ■■= ^TT-' + ^e>-p^^I)+p,,, (57) 

- p,,Z)+(i-^»^ + i-e>-p*^D+p,,) = ^vr«^Z)+p„, + i-eV-p^'^(Z?+p„fc)(Z)+ped). (58) 

Therefore, the trace of the equation (|54j) becomes 

11 1 

= 2A - n'^'D^pab + -e-^^/^.TTf, + -e->,f.p,rf7r'^%''^ + —Tr'^^D+pab 

= 2A + ie-^V'^^.VTb + ^"''( - D^pab + \e-'' PacPbdT^'"' + ^D+pab 

= 2A + ie-2>"^,7rb, (59) 
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where we used the equation (|1^. Thus A is given by 

A = -^e-^V^'^^TTeTTrf. (60) 

Thus we find that the equation (|54|) finahy becomes 

D-TT-' = -^e-^V-p'^^vr^vTrf + ^e'^'^ p-'p'^Tr.Tr^ - ie->,rf7r-7r^'^ 

+D+{^7r-' + ^eV-p^^lD+p,,)} + i-e>-pV^(^+/5ce)(Z?+Pd/). (61) 

One can show that this equation is the same as the equation H2.19p . To show this, let us multiply 
the equation (|HT|) by PamPbn-, and use the definitions of the conjugate momenta ^6.'6\ . • • •, (|3.5p . 
Then each term in (|61|) becomes as follows, 

(i) PamPbniD-TT"^) 

= PamPbnD-[ - ^e>'^V''(^+Ped) + /ieV"^p^^(D_p,d)} 

= -ie"(Z)_(7)(Z)+p„„) + ]^e" p'\D^Prr,c){D+Pnd) + ]^e" p'\D_Pnc){D+Pmd) 

-U^'iD^D+Pmn) + e''{D^h){D^p.mn) + he'' {D ^ a) (D ^ Pmn) 

-2h(fp'\D^Pmc){D^Pnd) + he^'iDlp^mn), (62) 

(ii) -e-^"^mT^n - -^e-^"" PmnP'^^T^c^d = ^e^'' PmcPndF+lF^^ - -e^" p^nPcdF+lF^t (63) 

(iii) J^e-''pcdTT'"'TT^'^ PamPbn 

= ^e''p''\D+Prac){D+Pnd) " ]^e'' p'\D+Pmc){D-Pnd) " ^e'^/>='^(D+p„e)(^-Pmd) 

+/ie>^'^(Z?_p^e)(I?-Pnd), (64) 

(iv) -Pa„^Pb„I?+{^v^'^'' + ^e>-p'«^(Z)+/5,,)} 

= -^Pam,06nI)+{e>'^V'''(i?-Pcd)} 

= -h''{D+a){D^p^n) + h"p'\D+p^,){D^Pr,d) + U''p'\D+Pnc){D-Pmd) 

-^e'^(Z)+i?_p„„), (65) 

(^) -^e>-p^'^/5^^(D+/5,e)(I?+Pd/)Pa™Pbn = -^e>'='^(D+p„,)(D+pw). (66) 

After a little algebra, we find that the equation (|6H) is identical to the equation (|2.19|) . Thus, 
assuming the boundary conditions Q, I have shown that the twelve Hamilton's equations of 
motion, (^ and 0, are just the first-order form of the six Einstein's equations (|2.15|) . H2.16() . 
H2.18() . and H2.19() . Therefore, the Hamilton's equations of motion, ^ and ©, together with 
the four divergence- type equations (|3.15jl . (|3.16jl . and (|3.17p . are completely equivalent to the full 
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Einstein's equations ^J^, ■■■, (|TT^ . 



[1] J.F. Plcbafiski, J. Math. Phys. 16 (1975) 2395. 

[2] Q-H. Park, Phys. Lett B 238 (1990) 287. 

[3] Q-H. Park, Phys. Lett B 257 (1991) 105. 

[4] Q-H. Park, Phys. Lett. B 269 (1991) 271. 

[5] L Bakas, Phys. Lett. B 228 (1989) 57. 

[6] E.G. Floratos, Phys. Lett. B 228 (1989) 335. 

[7] D.B. Fairhe, P. Fletcher, C.K. Zachos, J. Math. Phys. 31 (1990) 1088. 

[8] Y.M. Cho, Q.H. Park, K.S. Soh and J.H. Yoon, Phys. Lett. B 286 (1992) 251. 

[9] J.H. Yoon, Phys. Lett. B 308 (1993) 240. 

[10] J.H. Yoon, Phys. Lett. B 451 (1999) 296. 

[11] J.H. Yoon, Class. Quantum Grav. 16 (1999) 1863. 

[12] J.H. Yoon, J. Korean Phys. Soc. 34 (1999) 108. 

[13] J.H. Yoon, Phys. Lett. A 292 (2001) 166. 

[14] J.H. Yoon, (l+l)-dimensional formalism and quasi-local conservation equations, Advances in General 

Relativity and Cosmology, ed. G. Ferrarese (Pitagora Editrice Bologna 2002). 

[15] L.C. Evans, Partial Differential Equations, Graudate Studies in Mathematics 19 (American Mathemat- 
ical Society 1998). 

[16] J.D. Brown and J.W. York, Phys. Rev. D 47 (1993) 1407. 

[17] J.D. Brown, S. R. Lau, and J.W. York, Phys. Rev. D 59 (1999) 064028. 

[18] L.D. Landau and E.M. Lifshitz, The Classical Theory of Fields (Pergamon Press 1989). 

[19] C.R. Prior, Proc. R. Soc. Lond. A 354 (1977) 379. 

[20] R. Geroch, in Asymptotic Structure of Spacetime, ed. F.P. Esposito and L. Wittcn (Plenum Press, New 

York, 1977). 

[21] B.D. Bramson, Proc. R. Soc. Lond. A 364 (1978) 383. 

[22] A. Ashtekar and M. Streubel, Proc. R. Soc. Lond. A 376 (1981) 585. 

[23] T. Dray and M. Streubel, Class. Quantum Grav. 1 (1984) 15. 

[24] A. Ashtekar, L. Bombelli, and O. Reula, in Mechanics, Analysis and Geometry: 200 Years after La- 
grange, ed. M. Francaviglia (Elsevier 1991). 

[25] B. Carter, Comm. Math. Phys. 10 (1968) 280. 

[26] L. Hughston, R. Penrose, P. Sommers, and M. Walker, Comm. Math. Phys. 27 (1972) 303. 

[27] B. Carter and R.G. McLenaghan, Phys. Rev. D 19 (1979) 1093. 

[28] F. de Felice, J. Phys. A 13 (1980) 1701. 

[29] K. Kuchaf, Time and Interpretations of Quantum Gravity, Proceedings of the 4th Canadian Conference 

on General Relativity and Relativistic Astrophysics (World Scientific 1992). 

[30] J. Hartle, Astrophys. J. 157 (1970) 1395. 

[31] R. Wald, General Relativity (The University of Chicago Press 1984). 

[32] R. Penrose and W. Rindler, Spinors and Space-time, vol. 2 (Cambrige University Press 1986). 

[33] P.T. Chrusciel, J. Jezierski, and J. Kijowski, Hamiltonian Field Theory in The Radiating Regime 

(Springer- Verlag 2002). 

[34] Y.M. Cho, J. Math. Phys. 16 (1975) 2029. 



36 



[35] Y.M. Cho and P.G.O. Freund, Phys. Rev. D 12 (1975) 1711. 

[36] R.A. d'Inverno and J. Stachel, J. Math. Phys. 19 (1978) 2447. 

[37] R.A. d'Inverno and J. Smallwood, Phys. Rev. D 22 (1980) 1233. 

[38] E.T. Newman and T. Unti, J. Math. Phys. 3 (1962) 891. 

[39] E.T. Newman and K.P. Tod, Asymptotically Flat Space-Times in General Relativity and Gravitation, 

vol.2, ed. A. Held (Plenum Press, New York, 1980). 
[40] L. Tamburino and J. Winicour, Phys. Rev. 150 (1966) 1039. 
[41] J.H. Yoon, in preparation. 
[42] R. Abraham, J.E. Marsden, and T. Ratiu, Manifolds, Tensor Analysis, and Applications, Applied 

Mathematical Sciences 75 (Springer- Verlag 1983). 
[43] H. Bondi, M. van der Burg, and A. Metzner, Proc. R. Soc. Lond. A 269 (1962) 20. 
[44] J.D. Bekenstein, Phys. Rev. D 7 (1973) 949. 
[45] A.J. Chorin and J.E. Marsden, A Mathematical Introduction to Fluid Mechanics, Texts in Applied 

Mathematics 4 (Springer- Verlag 1992). 
[46] S. Frittelli, Well-posed first-order reduction of the characteristic problem of the linearized Einstein 

equations, ^gr-qc/6408035| 



